Coherent control of atomic tunneling 
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We study the tunneling of a two-level atom in a double well potential while the atom is coupled to 
a single electromagnetic held mode of a cavity. The coupling between internal and external degrees 
of freedom, due to the mechanical effect on the atom from photon emission into the cavity mode, 
can dramatically change the tunneling behavior. We predict that in general the tunneling process 
becomes quasiperiodic. In a certain regime of parameters a collapse and revival of the tunneling 
occurs. Accessing the internal degrees of freedom of the atom with a laser allows to coherently 
manipulate the atom position, and in particular to prepare the atom in one of the two wells. The 
effects described should be observable with atoms in an optical double well trap. 

PACS numbers: 73.40.Gk, 37.30.+i 



INTRODUCTION 



The tunneling effect is considered one of the hallmarks 
of quantum mechanical behavior. Historically, tunneling 
was first examined for single particles (e.g. a particles 
electrons in field emission [2 land later in mesoscopic cir- 
cuits J3|), for Cooper pairs J], and for molecular groups 
[S S0|- Recently the tunneling of atoms has attracted 
substantial attention 0, [^, [l(J, E[ • Dynamical (chaos as- 
sisted) tunneling of ultracold atoms between different is- 
lands of stability in phase space was analyzed in [l^, fl3j 
and has been observed experimentally [3, EH- Reso- 
nantly enhanced tunneling of atoms between wells of a 
tilted optical lattice has also been observed very recently 
[l6j |. fn all of these examples, the atoms have been con- 
sidered internally as inert, and only the center of mass 
coordinate of the atom was of interest. In [17I | it was 
shown that by taking into account the internal degrees 
of freedom of atoms, an atom/optical double well poten- 
tial could be created in which tunneling atoms see their 
internal and external states correlated (such an effect is 
also known from other contexts [HI). Mechanical effects 
of light in optical resonators were also investigated in 
[l9| . but no tunneling was considered. 



Here we show that the tunneling effect can be drasti- 
cally modified if an internal transition of the atom is cou- 
pled to a single electromagnetic mode in a cavity, such 
that photon emission is a reversible and coherent pro- 
cess. The resulting Rabi oscillations between states with 
the excitation in the atom and states with a photon in 
the cavity modulate the periodic tunneling motion. De- 
pending on the frequencies involved, a rich quasi-periodic 
behavior can result. If the cavity is fed with a coherent 
state, collapse and revival of the tunneling effect can oc- 
cur. Moreover, we show that one may profit from access 
to the internal degrees of freedom of the atom (e.g. with 
a laser) to control the atomic motion in the external po- 
tential. 





FIG. 1: (Color online) Two-level atom in a double well po- 
tential interacting with a standing wave inside a cavity. 



II. MODEL 

A. Derivation of the Hamiltonian 

Consider a trapped two-level atom (with levels \g), |e) 
of energy =fhu>o/2 respectively) interacting with a stand- 
ing wave (with wave number k and frequency uS) inside 
a cavity as illustrated in Fig. [TJ The atom is assumed to 
be bound in the y — z plane at the equilibrium position 
y — z — and to experience a symmetric double well po- 
tential V(x) along the x direction. We denote by A the 
tunnel splitting, i.e. the energy spacing between the two 
lowest energy states (the symmetric |— ) and antisym- 
metric |+) states) of this double well potential. Below 
we also allow the trapped atom to interact resonantly 
with an external laser. The Hamiltonian of this system 
is given by 



H = H A + H F + H AF , 



(1) 



where Ha 



H™ is the Hamiltonian of the trapped 



atom, Hp is the Hamiltonian of the free field and Haf 
is the interaction Hamiltonian describing the atom-field 
interaction. We have 



H F = 
H AF 
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where d denotes the atomic dipole, 

E = E u e (a + ah sin(fc(x — x )) (3) 



-^7, where 



is the electric field operator, with E w = 

eo is the permittivity of free space, V the electromag- 
netic mode volume, xq the abscissa at the left cavity 
mirror (x < 0), and e the electric field polarization vec- 
tor. We have introduced the operators af 1 (resp. cr° x ) 
for i = X,y,z as the Pauli spin operators in the basis 
{| e )i (resp. {|+),|— )})■ The operator x stands for 
the center-of-mass position of the atom, p x is the conju- 
gate momentum along the x axis, m denotes the atomic 
mass, and a (a^) the annihilation (creation) operator of 
the cavity radiation field. 

We adopt the two-level approximation which consists 
of taking into account only the two lowest motional en- 
ergy states. This requires the Rabi frequency \J Ag 2 + S 2 
(with S = u> — u>q the detuning between the cavity field 
and the atomic transition frequencies) to be much smaller 
than the frequency gap A between the upper motional 
states and the ground state doublet (see Fig. [J) . Within 
this approximation, Hamiltonian H A X becomes 



HI 



HA 



(4) 



and the position operator takes the form x = t^^x* with 
b/2 = (+\x\—). We can form states that are mainly con- 
centrated in the left/right wells, 



(5) 



\L) = (1+) - h»/V2, 
\R) = (]+) + \-))/V2. 



The average position of a particle localized in the right 
well is then given by b/2 (see Fig. [T]) and cr° x = |i£)(-R| — 
\L)(L\. The interaction Hamiltonian Haf can then be 
written 



Haf = — fig (a + a ) sin \ cos n erjf — cos x sin k o a *o™ 

with the atom-field coupling strength g 
— (e\d\g}.eE Ll j/h, and 



X = kxo, k = kb/2. 



(0) 



For long wavelengths (k < 1), or k = nix with inte- 
ger n, the left and right sites of the double well are in- 
distinguishable to the cavity photon and Haf reduces 
to Jaynes-Cummings Hamiltonian without rotating wave 
approximation (with a sine varying coupling constant), 
— hg sin x {a + a<)o*£: Note that n <C 1 would normally 
be identified with the Lamb-Dicke regime. Here the sit- 
uation is more subtle as the level spacings between the 
tunneling split ground state doublet and the next excited 
states can be very different such that the recoil energy 
Su'recoii satisfies A <C Lu IC con <C A. One may thus be in 



the Lamb-Dicke regime concerning transitions to higher 
vibrational states but have a significant mechanical ef- 
fect on the atomic tunneling. Furthermore, since there is 
only one photon mode, the recoil energy cannot vary con- 
tinuously and exciting higher vibrational levels requires 
^recoil close to a level spacing. Our numerical calculations 
show that even for k ~ 1 the two-level approximation can 
still work very well (see Fig. [|]). 

For (5, A <C u>, loo, a rotating wave approximation is 
justified, which consists in eliminating the energy non- 
conserving terms aa^cr™ and a^a^a™ with a™ = \e)(g\, 

a™ = and a c * = |+)(-|, cr ra = cr+ f . Within this 
approximation, the total Hamiltonian reads 



H 



hA 



?W0 in . t + 
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+%(aer" 1 + aV™) cos x sin k c^ x — sin x cos n l c 



(7) 



Thus, depending on the parameters x an d K i the cavity 
photon may induce internal transitions in the atom only 
(cos x sin k = 0), or induce transitions between internal 
and external states at the same time (cos x sin k ^ 0) 
even for a vanishing detuning (5 = to — = 0). This is 
in contrast to conventional sideband transitions of har- 
monically bound atoms or ions in the Lamb-Dicke regime 
which require an appropriate value of the detuning. For 
a fixed potential center (and thus fixed x)j K can be 
changed through a modulation of the well-to-well sep- 
aration b. We will neglect in the following the effects 
of decoherence, which means that not only g but also 
A should be much larger than the rate of spontaneous 
emission T, and the cavity decay rate K raT . 

We denote the global state of the atom-field system by 
\n,i,j) = |n)<g>|«)<8>|j) where \n) stands for the cavity field 
eigenstates, |i) £ {|— ),|+)} for the external motional 
states, and \j) £ {|g),l e )} for the internal states. The 
total excitation number N is given by a) a + crVV™. 



B. Energy levels 

The states |0, ±,g) are eigenstates of H with eigen- 
value (—ftwQ ± fiA)/2, i.e. these states remain uncou- 
pled and represent the two lowest energy states in the 
regime <5, A <C to, loq. It is straightforward to ver- 
ify that the Hamiltonian ([7|) only induces transitions 
between states with the same number of excitations 
N, {\N - l,+,e),\N,+,g),\N ~ l,-,e),\N,-,g)} = 
{|1), |2), |3), |4)}. It is therefore sufficient to solve the 
dynamics in this subspace. In doing so, we obtain the 
eigenvalues of H, 



Ml,, 



for p, p £ {±}, N = 1, 2, . . ., and with 



(8) 
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Q± = v / 2Ng 2 (l - cos(2k) cos(2x)) + S 2 + A 2 ± 2ft 2 , 

fl 2 = J 4Ng 2 cos 2 k sin 2 x(A 2 + 4Ng 2 sin 2 k cos 2 x) + <^ 2 A 2 • 



(9) 
(10) 



For a vanishing tunnel splitting (A = 0), ft± reduces to 
the maximum (minimum) of the two Rabi frequencies of 
the Jaynes-Cummings models in the right and left wells. 
For cosk = 1, the decoupling of external and internal 
degrees of freedom manifests itself also in the eigenvalues 
with ft± = W^Ng 2 sin 2 x + £ 2 ± A|. 



C. Evolution operator 

The whole dynamics of the system can be described 
by means of the evolution operator U(t) = e -' lHt / h with 
components Uj — (i\U(t)\j) = Uji, which can be calcu- 
lated exactly. In order to simplify the expressions, we 
restrict ourselves in the following to \ — — tt/4 — 2mr (in- 
teger n). We find, up to a an overall phase e -~ t ( N - 1 / 2 )^ t ; 

e/u = — ^ [iiS,si-^ {i + m(a - s)n 2 } 

m=± (ii) 

- vtft+ft_C M (<5A - fifl 2 )} , 



III. INTERNAL AND EXTERNAL DYNAMICS 

The reduced density matrix p cyi for the atomic center- 
of-mass motion alone follows from p = \?p(t)}(ip(t)\ by 
tracing out the field and internal degrees of freedom, 
where the total wave function at time t reads \ip(t)) = 
J2i i=i Uij(j\ip(0)} \i). The average position of the atom 
in the double well potential is then given by 



i\l Ng cos k 
V2A 



[pS^-^A 2 + 2Ng 2 sin 2 k 



+ fiil 2 ) + i/ift + ft_ AC P 



(12) 



Ui 3 = 



-iNg 2 sin(2/t) 
2A 



(13) 



with 



£ = A(S 2 + 2Ng 2 cos 2 k - 5A), 

a = ft, ft_ft 2 , 



(14) 
(15) 



and where all time dependence is in the coefficients 

C± = cos(ft±t/2), S ± = sin(ft ± i/2). (16) 

The remaining components can be deduced from the 
relations U 22 (S,A) = U 33 (-6,-A) = U U (S,-A) = 
U n (-S,A), U 24 (S,A) = U 13 (-S,A), Uu(S,A) = 
U 2 3(8,-A) = U 23 (-6,A), and U 34 (S,A) = U 12 (5,-A), 
valid for any where we have made explicit the depen- 
dence of the Uu on 8 and A. 



Tr ex (p e 



2 P ll) 



(17) 



with pll = {L\p cyi \L). Similarly, we obtain the reduced 
density matrix p ln for the internal atomic state by tracing 
out the field and external degrees of freedom, and the 
probability to find the atom in the excited state as p ee = 
(e|p in |e). 

In the following, we first focus on resonant atom-field 
interaction (uj — luq) before moving to the non-resonant 
case (a; ^ ujq). We distinguish three regimes according 
to the tunnel splitting compared to the Rabi frequency 
g : the small tunnel splitting regime (when A/g <C 1), 
the intermediate regime (when A/g ~ 1), and the large 
tunnel splitting regime (when A/g ^> 1). 



A. Resonant atom-field interaction 

For resonant atom-field interaction (5 = 0), the ex- 
pressions for Uj can be greatly simplified. If the system 
is initially prepared in the state \N — l,R,e) and for 
k = 7r/4, we have 



U 23 = WN 9^^ £ [^S.^SA + 2Ng 2 cos 2 k - ptf)} 



Pll 



A 2 + Ng 2 
with the tunnel frequency 

1 



sin 



fttun^ 



o t un = -(fi+ + fi_), 



and 



1 E^=± K ~ A 2 ) cos(lV) + 4A 2 cos ( 



(18) 



(19) 



Q-j--n_ , 



8(Ng 2 + A 2 ) 



(20) 

The atom position oscillates with a single frequency 
fttun given by Eq. (fT!)]) , whereas p ee evolves with three in 
general incommensurable frequencies ft + , ft_, and (ft+ — 
ft_)/2 giving rise to a quasi-periodic signal. 

For A/g < 1, Eq. HH]) leads to p LL ~ (up to order 
(A/g) 2 ), indicating that tunneling is suppressed. This is 



4 



already obvious from {7J , as the term responsible for tun- 
neling, {hA/2)cr c z x = (HA/2)(\R)(L\ + \L)(R\) becomes 
very small compared to the last term, diagonal in \R), \L) 
which leads to internal Rabi flopping. Note, however, 
that tunneling is suppressed on all time scales, even for 
t ^> 1/A, due to the reduced amplitude in Eq. (|18|) . very 
much in contrast to tunneling without internal degrees of 
freedom, where only the period of the tunneling motion, 
but not the amplitude is affected when A is reduced. For 
k approaching 7r, the situation changes because the term 
g cos x sin n cr® x of the interaction Hamiltonian inducing 
transitions between vibrational states becomes small in 
comparison with A thereby allowing tunneling again. 

Because internal and external degrees of freedom are 
coupled, the tunneling frequency (Eq. (|19j) ) depends on 
the number of photons inside the cavity. As an exam- 
ple, let us now consider A ~ g and a cavity field initially 
in a coherent state |a) = e^M J2"^=o \ n ) w i tn \ a \ 2 
equal to the mean photon number (n) . Figure [2] shows 
that the average position of the atom in the double well 
as a function of time for a coherent state exhibits col- 
lapses and revivals. The oscillation amplitude decreases 
with increasing mean photon number (n) and decreasing 
tunnel splitting A (see Eqs. ()18|9p ). Since the probabil- 
ity to find the atom in the excited state oscillates with 
three frequencies, no collapses and revivals are observed 
for p ee . 

The collapse time t c of the tunneling motion can be 
estimated from the condition Q {^tnn((n) + \J ( n )) — 

fitun((n) - \A™))) *c ~ 1 with Qtun(m) given by Eq. (JTSJ) 
for N = m + 1, which yields, for (n) 3> 1, 




FIG. 2: (Color online) Average position of the atom in the 
double well as a function of time for A/g — 2 (blue, top curve) 
and A/g = 5 (red), k = 7r/4 and a coherent state with a = 5. 
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(A/g) 2 + 3/4 
2(n) 



(21) 



The time interval between two following revivals, t r , 
follows from (f2tun((^)) — ^tun((") — 1)) t r = 2ir, and is 
given for (n) 3> 1 by 




O«n)"0 (22) 



For the parameters of Fig. [21 Eq. dHJ) yields gt r ~ 68.23 
for A/g — 2 and gt r ~ 86.70 for A/g = 5. Smaller revival 
times are possible for smaller values of (n), but in general 
the observation of revivals will be quite challenging, as 
they require A ~ g ^> K cav . 

For large tunnel splitting, A/g » 1, f2 tun = A + 
Ng 2 /(2A) + 0{(g/A) 3 ), and Eq. (JIgJ) reduces to p LL ~ 
sin 2 (At/2), which is identical to the tunneling of a parti- 
cle without internal structure. Equation (j2"D)) reduces to 
a Rabi oscillation p ee — cos 2 (\/Ngt/ 2). 



B. Non- resonant atom- field interaction 

For non- resonant atom-field interaction (5 ^ 0), 
intermediate tunnel splitting [see Fig. [3] for A = 5 - 



and 

-g], 



FIG. 3: (Color online) Average position of the atom in the 
double well as a function of time for A — 8 — g, k = 7t/4 and 
N — 1. The blue solid/red dashed curve corresponds to an 
excited atom initially located in the left/right well. 



(x(t)) involves in general the two non-commensurate fre- 
quencies f2+ and f2_ and varies therefore quasiperiod- 
ically as a function of time. Figure [3] also shows that 
an atom initially located in one of the two wells remains 
mostly confined to that well. 

For small tunnel splitting, A/g -C 1 and large detuning 
\5\/g > 1 (with A|<5|/g 2 ~ 1), the matrix elements of U 
simplify to 



Ul3 = 



u 33 



iNg 2 sin 2k 



sm 



S 2 A 2 + N 2 g 4 sin 2 (2k) 
Clt\ iSA 



ttt 



5 2 A 2 + N 2 g 4 sin 2 (2f. 



(23a) 
fit" 



2 

(23b) 



up to corrections of order 0(A/g) and a phase factor 

e i[(Ng 2 /5+S)-(2N~l)u J ]t/2 wMe ^ components JJ l2 and 

U23 are of order 0(A/g). In this situation, the system 
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FIG. 4: (Color online) Density matrix elements prr (top) 
and pee (bottom) as a function of the interaction time gt 
for an initially excited atom located in the right well and 
for the parameters A/g ~ 0.3336, S/g = 3, k = 7r/4, and 
N = 1. Numerical results from the propagation of the time 
dependent Schrodinger equation with Hamiltonian {1} and 
rotating wave approximation are represented by circles and 
analytical results by solid curves. The time propagation was 
done with (h = m = 1) g = 0.01 and the double well potential 
V(x) = 0.08a; 4 - x 2 yielding a tunnel splitting A ~ 0.003336 
and a ratio A/ A /ip + P ~ 44.4 » 1. 



oscillates only between the two states \N — 1, - 
| N — 1 , — , e) with a single frequency 



n 



S 2 A 2 +N 2 g A s i n 2(2 K ) 
6 



and 



(24) 



just as a three-level atom undergoing a Raman transition 
in the far detuned regime behaves as a two-level system. 

If the system is initially in the state \N — 1, — , e), we 
have from Eqs. (|23| 



(25) 



1 NS Asin(2«) r mA , 
PLL= 2- 20W t 1 - 008 ("*)]» 

and p ee = 1. For a detuning <5 = ±Ng 2 sin(2fc)/A, 



Pll 



1 1 

2 T 4 



1 



(\/2Ai) 



(26) 



This regime may be suitable for coherently manipu- 
lating the atom position through access to its internal 
degrees of freedom with a laser. Coherent manipulation 
of the position of neutral atoms has been proposed and 
demonstrated before, see e.g. [U [53, HH . In these exam- 
ples, the manipulation is done by modifying the external 
potential. The mechanism we propose here is very dif- 
ferent, as the potential remains totally unchanged, and 
only internal transitions and the tunneling effect are used 
to move the atom in a controlled way. As an example, 
we show how the atom can be prepared in the left well 
starting from the ground state |0, —,g) for 5 — —g 2 /A. 

We first apply a 7r-pulse with an external laser resonant 
with the atomic transition. By using a laser with a wave 



vector perpendicular to the Ox-direction, only the atomic 
internal degree of freedom is affected, resulting in the 
transition \0,—,g) — > — i\0, — , e). We assumed that the 
laser Rabi frequency Qr is much larger than the tunnel 
frequency A. 

Now we use the coupling between the internal and ex- 
ternal degrees of freedom to create a superposition of 
the |0,±,e) states, and then apply a second resonant ir- 
pulse to get back to the uncoupled states \0,±,g). For 
A/g -C 1, S = —g 2 /A and K = 7r/4, the initial state 
transforms according to 



\o,-,g) 



O nt— 7r 



|0,-,e) p \0,L,e) — — \0,L,g) 

At=7r/v2 S2 R t=7T 

(27) 

up to a physically irrelevant phase. Other coherent su- 
perpositions of |0, +,<?) and |0, — ,g) can be obtained by 
choosing appropriate interaction times. 

In order to verify that the two-level approximation for 
the external motion used in the derivation of the Hamilto- 
nian is a good approximation, we have numerically solved 
the time dependent Schrodinger equation with Hamil- 
tonian ([1]) and rotating wave approximation but with 
the exact external potential V{x) (i.e. with a large num- 
ber of vibrational states). Figure |4] shows that provided 
A 3> y/ Ag 2 + S 2 as stated before, to take only the two 
lowest vibrational states into account is indeed a good 
approximation. 

We finally comment on possible experimental realiza- 
tions of our model. Double well potentials with tunable 
well-to-well separation have been demonstrated with op- 
tical dipole traps e.g. in [2l|, |2J], and on atom chips 
e.g. in [2^, [26[. For our model, the double well poten- 
tial has to be realized inside the cavity. Optical trapping 
and even cooling of atoms close to their ground state 
inside a cavity has been achieved in several groups by 
now [13, [H, H^, [10|, but up to our knowledge double 
well potentials have not been realized in a cavity so far. 
However, some of the cavities developed have a very long 
lateral opening (up to 222 pm [3l| ) and should allow more 
complicated trapping potentials (optical lattices inter- 
secting a cavity have been realized in Chapman's group 
[3l|). We remark that it is not essential for our model 
that the double well potential be aligned with the cavity 
axes. Any other orientation is possible, and only leads to 
modified coefficients cos x sin k and sin \ cos k. 

At certain "magical wavelengths" , Cs, Yb, Sr, Mg, and 
Ca atoms in optical traps experience the same potential 
for ground and excited internal states coupled by a dipole 
transition [53, H2, [H, H3| ■ In a symmetric potential V(z) 
the tunneling frequency A is given in WKB approxi- 
mation by A ~ uj osc exp(— l/h ^2m(Eo — V(z)) dz) 
where Eq is the ground state energy, uj osc the single well 
harmonic oscillation frequency, and z = ±a are the corre- 
sponding classical turning points delimiting the range of 
the barrier. The exponential factor can approach unity 
for a barrier that is only slightly higher than the ground 
state energy Eq, in which case cooling to temperatures 
fcgT < HA should be possible with state of the art 
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techniques [27j. In [27] a trap depth V /h = 47 MHz 
was achieved inside a cavity with 1.2 mW laser power. 
In any case, the trap frequency and thus the tunneling 
splitting are determined by the laser power and the fo- 
cussing (or the wavelength for optical lattices), and can 
therefore be controlled independently of T, K C av, such 
that there should be no fundamental problem achieving 
A 3> r, K cav . The detection of the tunneling motion 
should be possible by optical imaging, i.e. diffusion of 
laser light from another transition in the optical regime 
with smaller wavelength than the well separation. Al- 
ternatively, one might monitor the transmission through 
the cavity in the case that it differs for the two locations 



of the wells [35| . Another possibility might be using the 
atomic spin as a position meter [l7j |. 
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